Hyperfine interaction in a quantum dot: Non-Markovian electron spin dynamics 
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We have performed a systematic calculation for the non-Markovian dynamics of a localized elec- 
tron spin interacting with an environment of nuclear spins via the Fermi contact hyperfine inter- 
action. This work applies to an electron in the s-type orbital ground state of a quantum dot or 
bound to a donor impurity, and is valid for arbitrary polarization p of the nuclear spin system, 
and arbitrary nuclear spin / in high magnetic fields. In the limit of p = 1 and J = |, the Born 
approximation of our perturbative theory recovers the exact electron spin dynamics. We have found 
the form of the generalized master equation (GME) for the longitudinal and transverse components 
of the electron spin to all orders in the electron spin-nuclear spin flip-flop terms. Our perturbative 
expansion is regular, unlike standard time-dependent perturbation theory, and can be carried-out 
to higher orders. We show this explicitly with a fourth-order calculation of the longitudinal spin 
dynamics. In zero magnetic field, the fraction of the electron spin that decays is bounded by the 
smallness parameter 8 = l/p 2 N, where N is the number of nuclear spins within the extent of the 
electron wave function. However, the form of the decay can only be determined in a high magnetic 
field, much larger than the maximum Overhauser field. In general the electron spin shows rich 
dynamics, described by a sum of contributions with non-exponential decay, exponential decay, and 
undamped oscillations. There is an abrupt crossover in the electron spin asymptotics at a critical 
dimensionality and shape of the electron envelope wave function. We propose a scheme that could 
be used to measure the non-Markovian dynamics using a standard spin-echo technique, even when 
the fraction that undergoes non-Markovian dynamics is small. 
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I. INTRODUCTION 

Prospects for the development of new spintronic 
devices^ and the controlled manipulation of electron or 
nuclear spins for quantum information processing^ have 
sparked substantial research efforts in recent years. One 
of the major obstacles to achieving these goals is deco- 
herence due to the influence of an uncontrollable environ- 
ment. For quantum computing tasks, the strict require- 
ments for error correction^ put strong limits on the degree 
of decoherence allowed in such devices. From this point 
of view, single-electron semiconductor quantum dots rep- 
resent good candidates for spin-based information pro- 
cessing since they show particularly long longitudinal re- 
laxation times, T\ = 1 ms. 4 In GaAs quantum wells, the 
transverse dephasing time T 2 * f° r an ensemble of electron 
spins, which typically provides a lower bound for the in- 
trinsic decoherence time T2 of an isolated spin, has been 
measured to be in excess of 100 ns£ 

Possible sources of decoherence for a single electron 
spin confined to a quantum dot are spin-orbit coupling 
and the contact hyperfine interaction with the surround- 
ing nuclear spins^ The relaxation rate due to spin-orbit 
coupling jr is suppressed for localized electrons at low 
temperatures 7 ^ and recent work has shown that T2, due 
to spin-orbit coupling, can be as long as T\ under re- 
alistic conditions^ However, since spin-carrying isotopes 
are common in the semiconductor industry, the contact 
hyperfine interaction (in contrast to the spin-orbit inter- 
action) is likely an unavoidable source of decoherence, 
which does not vanish with decreasing temperature or 



carefully chosen quantum dot geometry^ 

In the last few years, a great deal of effort has been 
focused on a theoretical description of interesting ef- 
fects arising from the contact hyperfine interaction for 
a localized electron, 6 ! 11 ! 12 ! 13 ! 14 ! 15 ! 16 ! 17 ! 18 ! 19 ! 20 ! 21 ! 22 ! 23 ! 24 ! 25 
The predicted effects include a dramatic variation of T\ 
with gate voltage in a quantum dot near the Coulomb 
blockade peaks or valleys^ all-optical polarization of 
the nuclear spins^ use of the nuclear spin system as 
a quantum memory 1 ^ and several potential spin re- 
laxation and decoherence mechanisms ,JiI2iiiL2L2I This 
theoretical work is spurred-on by intriguing experiments 
that show localized electrical detection of spin resonance 
phenomena^ nuclear spin polarization near quantum 
point contacts^ 7 , gate-controlled transfer of polarization 
between electrons and nuclei^ nuclear spin polarization 
and manipulation due to optical pumping in GaAs quan- 
tum wells^ 9 . and voltage-controlled nuclear spin polar- 
ization in a field-effect transistor^ In addition, recent 
experiments have shown hyperfine induced oscillations 
in transport current through a double quantum dotj 3 ^ 
and long T2 times for electrons trapped at shallow donor 
impurities in isotopically purified 28 Si:Pi^ Our system of 
interest in this paper is an electron confined to a single 
GaAs quantum dot, but this work applies quite generally 
to other systems, such as electrons trapped at shallow 
donor impurities in Si:P^S 

In this paper, we investigate electron spin dynamics 
at times shorter than the nuclear dipole-dipole correla- 
tion time Tdd (Tdd ~ 10~ 4 s in GaAs is given directly 
by the inverse width of the nuclear magnetic resonance 
(NMR) line^ 3 .). At these time scales, the relevant Hamil- 
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tonian for a description of the electron and nuclear spin 
dynamics is that for the Fermi contact hyperfine inter- 
action (see Eq. 0), below). Dynamics under the ac- 
tion of this Hamiltonian may be of fundamental interest, 
since in zero magnetic field, Eq. corresponds to the 
well-known integrable Gaudin magnet, which is soluble 
via Bethe ansatziifli^i Though the Hamiltonian appears 
simple, a detailed microscopic description for the dynam- 
ics of a spin coupled to a spin environment remains an 
open quest ion iSiSi A degree of success has been achieved 
some time ago in bulk systems through the development 
of phenomenological models^! These models invoke cer- 
tain approximations, namely, assumptions of Markovian 
dynamics and ensemble averaging. Care should therefore 
be taken in applying the same models to the problem of 
single-spin dccohercncc for an electron spin strongly cou- 
pled to a nuclear spin environment, where they may not 
applyi n i 12 

For nuclear spin I = », an exact solution for the 
electron spin dynamics has been found in the special 
case of a fully polarized initial state of the nuclear spin 
systemiii*i2 This solution shows that the electron spin 
only decays by a fraction oc of its initial value, where 
N is the number of nuclear spins within the extent of the 
electron wave function. The decaying fraction was shown 
to have a non-exponential tail for long times, which sug- 
gests non-Markovian (history dependent) behavior. For 
an initial nuclear spin configuration that is not fully po- 
larized, no exact solution is available and standard time- 
dependent perturbation theory failsiii Subsequent ex- 
act diagonalization studies on small spin systems^ have 
shown that the electron spin dynamics are highly de- 
pendent on the type of initial nuclear spin configuration, 
and the dynamics of a randomly correlated initial nu- 
clear spin configuration arc reproduced by an ensemble 
average over direct-product initial states. The unusual 
(non-exponential) form of decay, and the fraction of the 
electron spin that undergoes decay may be of interest 
in quantum error correction (QEC) since QEC schemes 
typically assume exponential decay to zero. 

In this paper we formulate a systematic perturbative 
theory of electron spin dynamics under the action of the 
Fermi contact hyperfine interaction. This theory is valid 
for arbitrary nuclear spin polarization and arbitrary nu- 
clear spin / in high magnetic fields. For nuclear spin 
I = \ and a fully polarized nuclear spin system, we re- 
cover the exact solution for the electron spin dynamics 
within the Born approximation of our perturbative the- 
ory. Our approach follows a method recently applied 
to the spin-boson models This method does not suf- 
fer from unbounded secular terms that occur in standard 
perturbation theory^ and does not involve Markovian 
approximations. 

This paper is organized as follows. In Section [H] we 
review the model Hamiltonian and address the question 
of realistic initial conditions. In Section ITTll we derive the 
form of the exact generalized master equation (GME) 
for the electron spin dynamics. In Section IIVI we con- 



sider the leading-order electron spin dynamics in high 
magnetic fields. In Section[V]we proceed to calculate the 
complete non-Markovian dynamics within the Born ap- 
proximation. We describe a procedure that could be used 
to measure the non-Markovian dynamics in Section IVII 
In Section IVlII we show that our method can be extended 
to higher orders without the problems of standard per- 
turbation theory by explicitly calculating the corrections 
to the longitudinal spin self-energy at fourth order in the 
nuclear spin-electron spin flip-flop terms. We conclude 
in Section lYlIII with a summary of the results. Technical 
details are deferred to Appendices IATIeI 

II. MODEL 

A. Hamiltonian 

We consider a localized electron spin interacting with 
iVtot nuclear spins via the Fermi contact hyperfine inter- 
action. The Hamiltonian for this system is 

H = bS z + e nz I z +h-S, (1) 

where S = (S x , S y , S z ) is the electron spin operator, b = 
g* [ibB z (e nz = gifJ,NB z ) is the electron (nuclear) Zeeman 
splitting in a magnetic field B z , with effective g- factor g* 
(gi) for the electron (nuclei) and Bohr (nuclear) magne- 
ton li b (fj, N ). Further, h = {h x ,h y ,h z ) = X^=o _1 A k^k 
gives the (quantum) field generated by an environment 
of nuclear spins. Ifc = is the nuclear spin op- 

erator at lattice site k and is the associated hyperfine 
coupling constant. I z = ^ k If. is the total z-componcnt 
of nuclear spin. 

The nuclear Zeeman term can be formally eliminated 
from the Hamiltonian H (Eq. QJ) by transforming to a 
rotating reference frame. The z-component of total angu- 
lar momentum is J z = S z + I z . Adding and subtracting 
tnzJz gives Ti = TV + e nz J z . The Hamiltonian in the 
rotating frame, TV , is then 

Ti = TV + H' v , (2) 
H' = b'S z + h z S z , (3) 

H'y = ±(h + S^+h^S + ), (4) 

where b' = b—e nz and we have introduced h± = h x ±ih y . 
The usual Heisenberg-picture operators in the rotating 
frame are S' x (t) = e mt S x e- i ' H ' t , X = z,+, S± = S x ± 
iS y . Noting that [J Z ,TV\ = 0, we find they are related to 
the operators Sx(t) = e '™5xe~ in the rest frame by 

S' z (t) = S z (t) (5) 
S' + (t) = e-^S+it). (6) 

In the following, (S' z ) t and (S' + ) t will be evaluated in the 
rotating frame, but we omit primes on all expectation 
values. 
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Figure 1: Schematic of the square modulus of the electron 
envelope wave function |?/>(j")| 2 and nuclear spins (arrows), k 
is the nuclear site index, N is the number of nuclear spins 
within radius r = lo, and iVtot is the total number of nuclear 
spins in the system. 

The hyperfine coupling constants A k are given byiS 

A k = Av \^(r k )\ 2 . (7) 

Here, vq is the volume of a crystal unit cell containing 
one nuclear spin, tp( r ) is the electron envelope wave func- 
tion, and A is the strength of the hyperfine coupling. In 
GaAs, all naturally occurring isotopes carry spin I = ^. 
In bulk GaAs, A has been estimated^ to be A = 90 p,eV 
( ,-^4 — = 3.5 T) . This estimate is based on an average 

over the hyperfine coupling constants for the three nu- 
clear isotopes 69 Ga, 71 Ga, and 75 As, weighted by their 
relative abundances. Natural silicon contains 4.7% 29 Si, 
which carries I = \, and 95% 28 Si, with 1 = 0. An 
electron bound to a phosphorus donor impurity in natu- 
ral Si:P interacts with N as 10 2 surrounding 29 Si nuclear 
spins, in which case the hyperfine coupling constant is on 
the order of A m 0.1 fieVML We consider a localized elec- 
tron in its orbital ground state, described by an isotropic 
envelope wave function of the form 



ip{rk) = ip(0)exp 







H 





(8) 



When m — 2, ip(r) is a Gaussian with Bohr radius Zrj, and 
for m = 1, ip(r) corresponds to a hydrogen- like s-state 
with Bohr radius ao = 2lo- N tot nuclear spins are in the 
system, but the effective number N of spins interacting 
appreciably with the electron is smaller (see Fig. 2J. N 
is defined as the number of nuclear spins within radius lo 
of the origin and the integer index k gives the number of 

spins within radius r k . In d dimensions, (j^) = It 

= 1, 
0). 



is convenient to work in energy units such that 
where Aq is the coupling constant at the origin (tq - 
In these units A k takes the simple form 



A k = 2 exp 



(9) 



B. Initial conditions 

1. Sudden approximation 

The electron spin and nuclear system are decoupled 
for times t < 0, and prepared independently in states 
described by the density operators ps{0) and pi(0), re- 
spectively. At t = 0, the electron and nuclear spin system 
are brought into contact "instantaneously", i.e., the elec- 
tron spin and nuclear system are brought into contact 
over a switching time scale t sw ^ which is sufficiently 
small-see Eq. I)llfl . below. The state of the entire sys- 
tem, described by the total density operator pit) is then 
continuous at t = 0, and is given by 



p(0-)=p(0 + ) = Ps (0)® Pl (0). 



(10) 



The evolution of the density operator p(t) for t > is gov- 
erned by the Hamiltonian TC for an electron spin coupled 
to an environment of nuclear spins. Since the largest en- 
ergy scale in this problem is given by \b' + A\. in general 
the condition 



2irh 



A\ 



(11) 



should be satisfied for the sudden approximation (Eq. 
|[TDJ) to be valid. In bulk GaAs, ^ ~ 50 ps and for an 
electron bound to a phosphorus donor in natural silicon, 
^ ~ 10 ns. 



2. Dependence on the nuclear state: zeroth order dynamics 

Evolution of the electron spin for different initial nu- 
clear co nfig urations has been addressed previously^Sii^ 
In Ref. UJ it was found, through numerical study, that 
the dynamics of the electron spin were highly dependent 
on the initial state of the nuclear system. The goal of 
this section is to shed more light on the role of the initial 
nuclear configuration by evaluating the much simpler ze- 
roth order dynamics, i.e., the electron spin evolution is 
evaluated under Ti' = H.' alone, neglecting the flip-flop 
terms TL' V . 

Since [H' ,S Z ] = 0, (S z ) t is constant. However, 
S±] 7^ 0, so the transverse components, (S+) t = 
f i(Sy} t , will have a nontrivial time depen- 
We evaluate the expectation value (S+) t = 

m o t S' + e J ' H o t p(0)| (setting h = 1), with the ini- 



K 

dence 
Tr < e 



tial state given in Eq. (|10|) . After performing a partial 
trace over the electron spin Hilbcrt space, we obtain an 
expression in terms of the initial nuclear spin state: 

(S + ) t = (S + ) o Tr I {e^' +h ^p I (0)}, (12) 

where Tij is a partial trace over the nuclear spin space 
alone. For simplicity, here we consider I = ^, and the 
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coupling constants are taken to be uniform. After en- 
forcing the normalization ^ fc = 2N in units where 
4 s1 = = 1) the hyperfine coupling constants are 



At 



2, jfc = 0,l---JV-l 
0, k > N 



(13) 



The zcroth-order electron spin dynamics can now be eval- 
uated exactly for three types of initial nuclear spin con- 
figuration: 



N 



P(N r ,N,ft)\N r )(Ni\ 

N r =0 

\n) („|. 



(14) 
(15) 

(16) 



(i) 



Pj ' is a pure state, where \tpi(0)} = 

YltoiVftlW+e^VT^hlik)) is chosen to 
render the z-component of nuclear spin translationally 
invariant: (V>/(0)| I z k |^/(0)) = ±(2/ T - 1) = § , and 
p = 2/j — 1 is the polarization of the nuclear spin 
system. <pk is an arbitrary site-dependent phase factor. 

P(x;n,f) = ( ™ ^ f x {\ — f) n ~ x is a binomial distribu- 
tion, and | TV-j- ) is a product state of the form ||TI " " " ) 

with iVf spins up and N — iV-j- spins down, (0) then 
corresponds to a mixed state; this is an ensemble of 
product states where the N spins in each product state 
are selected from a bath of polarization p = 2/f — 1. 

p^f \ like pj , is a pure state, but for this state |n) is 
chosen to be an eigenstatc of h z with eigenvalue pN 
(corresponding to a nuclear system with polarization p) : 
h z \n) = pN \n). We insert the initial nuclear spin states 
p^j\o) into (|T2*|) to obtain the associated time evolution 



(S. 



+/t 



AT 



7V T =0 



(S 



./! 3) 



(^+)o 



J(b'+pN)t 



(18) 



M(iVj) = 2Aj — A is the nuclear magnetization on a dot 
with Nf nuclear spins up. 

The similarity in dynamics between randomly corre- 
lated (entangled) pure states and mixed states has been 
demonstrated for evolution under the full Hamiltonian 
(H' = H' + Hy) via exact diagonalizations of small 
(Atot ~ 19) spin systems^ Here, the zeroth order elec- 
tron spin dynamics are identical for the pure state p^f 1 (0) 

(2) 

and the mixed state p\ (0) even when the initial pure 
state \ipj(0)) is a direct product. Direct application of 
the central limit theorem gives a Gaussian decay for large 



N: 



(S. 



(S 



r+i(b'+ P N)t 



+ /0 ( 



(19) 

Returning to dimension-full units (c.f. Table [I] below), 
the time scale for this decay is given by t c = 2 -^-p£ c ~ 
5 ns for a GaAs quantum dot with p 2 <C 1 containing 
N = 10 5 nuclei and t c f=a 100 ns for an electron trapped 
at a shallow donor impurity in Si : P, with N = 10 2 . 
For an ensemble of nuclear spin states, Gaussian decay 
with the time scale r c has been found previouslyii^iSi 2 ^ 
Gaussian decay for a Hamiltonian with an Ising coupling 
of electron and nuclear spins has been demonstrated^ for 
a more general class of pure initial states and for coupling 
constants A k that may vary from site-to-site. 

f 1 2) 

For the initial states p\ ' (0), precise control over the 
nuclear spin polarization between measurements or a 
spin-echo technique would be needed to reduce or elim- 
inate the rapid decay described by (|19|) . However, the 
quantum superposition of h z eigenstates can be removed, 

in principle, from the pure state Pj (0) by performing 
a strong (von Neumann) measurement on the nuclear 
Overhauscr field pNm& After the nuclear system is pre- 
pared in an /i^-cigenstate, to zeroth order the electron 
spin dynamics will be given by (S+) t , i.e., a simple pre- 
cession about the z-axis with no decay. 

When higher-order corrections are taken into account, 
and the coupling constants A k are allowed to vary from 
site-to-site, even an initial /i z -eigenstate can lead to decay 
of the electron spin. This has been showniiiiSi in an exact 
solution for the specific case of a fully-polarized system 
of nuclear spins- ^ and by exact diagonalization on small 
systems^ The goal of the present work is to perform 
an analytical calculation with a larger range of validity 
(a large system of nuclear spins with arbitrary polariza- 
tion and arbitrary nuclear spin / in a sufficiently strong 
magnetic field) that recovers previous exact results in the 
relevant limiting cases. In the rest of this paper, the ef- 
fect of higher- (beyond zeroth-) order corrections will be 
considered for a nuclear spin system prepared in an ar- 
bitrary h z eigenstatc: pi(0) = p\ (0), as given in Eq. 
(|16fl . Specifically, the initial state of the nuclear system 
\n) can be written as an arbitrary linear combination of 
g n degenerate product states: 



JV t , 



(20) 



where |J, mj) is an eigenstate of the operator If with 
eigenvalue m, and h z \rij) = [ft«] nn \nj) for all j, where we 
write the matrix elements of any operator O as (i \ O \ j) = 

PL- 
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III. GENERALIZED MASTER EQUATION 

To evaluate the dynamics of the reduced (electron spin) 
density operator, we introduce a projection supcropcr- 
ator P, defined by its action on an arbitrary opera- 
tor O: PO = pi(0)TrjO. P is chosen to preserve all 
electron spin expectation values: (5 l g) t = Tr Spp(t) = 
Tr SpPp(t), (3 = x,y,z, and satisfies P 2 = P. For fac- 
torized initial conditions (Eq. JHHO, Pp{0) = p(0), 
which is a sufficient condition to rewrite the von Neu- 
mann equation p(t) = — i[H',p(t)] in the form of the 
exact Nakajima-Zwanzig generalized master equation 
(GME)22: 



Pp(t) = -iPLPp(t)-i / dt'E(t-t')p(t'), (21) 

Jo 

E(i) = -iPLQe- lLQt QLP, (22) 

where E(i) is the self-energy superoperator and Q = 1— P 
is the complement of P (1 is the identity superoperator). 
L = Lq + Ly is the full Liouvillian, where L a (a = 
V,0) is defined by L a O = [H' a ,0], When the initial 
nuclear state is of the form pi(0) = \n) (n|, where \n) is 
an arbitrary eigenstate of h z , as in Eq. I|2U|) , P obeys the 
useful identities 



Details of the expansion (Eq. (EHJ) arc given in 
Appendix ^ It is most convenient to evaluate the 
inhomogeneous term N z (t) and the memory kernels 
E zz (i), E++(t) in terms of their Laplace transforms: 
f(s) = J °° dte- st f(t), Re[s] > 0. N z (s) and E zz (s) are 
given in terms of matrix elements of the reduced self- 
energy by 



N z (s) = -_(E TT ( S )+E n ( S )), (31) 
E zz ( s ) = E TT (s)-S u ( S ). (32) 

Explicit expressions for the matrix elements 
E ++ (s), S||(s), and Ef^(s) are given in Appendix 
IaI We find that the self-energy at (2£;) th order is 
suppressed by the factor A fe , where 



N 



(33) 



The parameter A and some other commonly used sym- 
bols are given in dimensionless and dimension-full units 
in Tabic U below. For high magnetic fields \b'\ 3> N 



(\B Z \ 



jr\ <C 1, and the ex- 



PL V P = 0, 
PL P = L P. 



(23) 
(24) 

We apply Eqs. (|23|) and Q24JI. and perform a trace on 
1 12 11 1 over the nuclear spins to obtain 



we have |A| 

pansion is well-controlled. The non-perturbative regime 
is given by |A| > 1, and the pcrturbativc regime by 
| A | < 1. Thus, a perturbative expansion is possible when 
the electron Zeeman energy produced by the magnetic 
and/or Ovcrhauser field (provided by N nuclear spins) is 
larger than the single maximum hypcrfinc coupling con- 
stant A. In the rest of this section we apply the Born 



p S (t) = -il% Ps (t)-i 



f (2) 

/ dt'Es(i — t')ps{t')i (25) approximation Es ~ Eg to the reduced self-energy, and 
Jo perform the continuum limit for a large uniformly polar- 



£s(t) 



iTr I Le- lQLt L vPl (0), 



where U l O = [S z uj n ,0] and u n = b' + [h z ] nn . E<j(i) is 
the reduced self-energy superoperator. ps(t) = Ti jp(t) = 
9C0 + (S x ) t a x + (S y ) t <J y + (Sy) . <j y is the reduced electron 
spin density operator, where erg, /3 = x,y, z, are the usual 
Pauli matrices and ero is the 2x2 identity. 
We iterate the Schwingcr-Dyson identity^ 



perform the continuum limit for a large uniformly polar- 
(26) ized nuclear spin system. Later, we also consider higher 
orders. 



-iQ(L +L v )t 



e -' lQLot -i 



di . /e -iQL (t-t') QLve -iQLt' 



A. Born approximation 

In Born approximation, the memory kernels 
E zz (i), S + _)_(t) and inhomogeneous term N z (t) 
in l|29|) and l|3(J|) are replaced by the forms ob- 
tained from the lowest-order self-energy, i.e., 



on (|26|l to generate a systematic expansion of the reduced 
self-energy in terms of the perturbation Liouvillian Ly: 



(27) N z (t) -> 7vi 2) (i), E„(i) 



£&>(i), E ++ (t) 



(2), 



.(2) 



,(2) 
J + + 



E 5 (t) = 4 a) (t)+E^(t) + 



(4)/ 



(28) 



where the superscript indicates the number of occur- 
rences of Ly. Quite remarkably, to all orders in Ly, 
the equations for the longitudinal ((S z ) t ) and transverse 



(*)■ 

In Laplace space, £l?(s), E^| (s), and £++(s) are given 
for an arbitrary initial h z eigenstate |n) (see Eq. (I2UI0 
in Appendix El Eqs. ljA20)l . l|X2"T]l . and (JX22J). Inserting 
an initial state \n) for a large nuclear spin system with 
uniform polarization gives (see Appendix |B"|) : 



{(S. 



+/t 



x/ t 



(S y ) t ) electron spin components are 



decoupled and take the form: 

(S z ) t = N z (t)-i f dt'H zz {t-t l ){S z ) t 
Jo 

ft 



E t 



(29) 



,(2), 



(S+)t = iuJ n (S+) t - 



i 1 dt'Y; ++ (t-t')(S + ) t , .(30) 




I±(s) 



-iNc+ [I+(s - iw n ) + I-{s + iu n )] , (34) 
iNc_ [I_(s - iuj n ) + I+(s + iu> n )] , (35) 
-iN[c-I+(s)+c+I-(a)], (36) 

4A^7?%- (3?) 
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In the above, the coefficients 



c± =/(/+!)- «m(m±l))) 



(38) 



have been introduced, where ((F(m))) = 
Sm=-/ Pi{m)F{m) for an arbitrary function F(m). 
Pi(m) is the probability of finding a nuclear spin I 
with ^-projection m. The polarization p of the initial 
nuclear state is defined through the relation ((to)) = pi. 
Without loss of generality, in the rest of this paper 
p > 0, but b' may take on positive or negative values. 
Assuming a uniform polarization in the nuclear spin 
system, we can evaluate the nuclear Overhauser field in 
terms of the initial polarization: 



[h,] m = Y l M{^))=pIA, 



(39) 



where we have used J2i^i = A. 

The continuum limit is performed by taking iVtot —* 
oo, while TV ^> 1 is kept constant. For times t <C VN, 
this allows the replacement of sums by integrals ^ fc — > 
J °° dk, with small corrections (see Appendix . We 
insert the coupling constants Ak from Eq. ^ into Eq. 
(|37|l . perform the continuum limit and make the change 
of variables x = 4r- to obtain 



m J 



xlln: 



dx- 

o s =F ix 



v= 1. 

TO 



(40) 



We use the relation I± (t = 0) = lim s 
the initial amplitude 



J = J ± (i = 0) = ^-(i 
to \ 2 



s/±(s) to obtain 



TO 



(41) 



for an arbitrary ratio — . For parabolic confinement in 
two dimensions, m = d = 2. The integral in (|40|) can 
then be performed easily, which yields 



I±{s) = s[log(s=F«) -log(s)] ±i (m = d = 2). (42) 



In dimcnsionlcss units 



1, we find A = J2kAk 



J dkAk, with the coupling constants Ak given in Eq. 10: 



A = A N— r — = 2N—T — 



(43) 



IV. HIGH FIELD SOLUTION 

In the next section, we will obtain a complete solu- 
tion to the GME within the Born approximation. This 
complete solution will exhibit non-perturbative features 
(which can not be obtained from standard perturbation 
theory), in the weakly pcrturbative regime for the self- 
energy, which we define by |A| < 1. Here, we find the 



leading behavior in the strongly pcrturbative (high mag- 
netic field) limit, defined by |A| <C 1, or equivalently, 
1 6' | 3> N. We do this in two ways. First, we ap- 
ply standard perturbation theory, where we encounter 
known difficultiesii (secular terms that grow unbounded 
in time). Second, we extract the leading-order spin dy- 
namics from the non-Markovian remainder term in a 
Born-Markov approximation performed directly on the 
GME. We find that the secular terms are absent from 
the GME solution. We then give a brief description of 
the dependence of the spin decay on the form and dimen- 
sionality of the electron envelope wave function. 



A. Perturbation theory 

Applying standard time-dependent perturbation the- 
ory (see Appendix [DJ) to lowest (second) order in H' v , 
performing the continuum limit, and expanding the re- 
sult to leading order in -j- , we find 



(S+) t = 

(Sz) t = (^)oc + 



where 



,_osc 



(t) 



cc {t) 

h (*) 



cr soc 



[l-SI (c + + c-)} (S+) . 
5[C+I-(t)+Ctl+(t)], 
iAI {c+ + c-)(S+) t, 



(44) 
(45) 



(46) 
(47) 
(48) 



and 



(^)oo = [1 - 2 «o (c+ + c_)] (S Z ) Q + 2pUI , (49) 
af c {t) = 2<5Rc [e"^"* (C z + I_{t) + C* /+(*))] .(50) 

We have introduced the smallness parameter S = and 
the coefficients 



\ c± (5+) , X = +. 



(S z ) t is the sum of a constant contribution (Sz)^ and 
a contribution that decays to zero <xf ec (i) with initial 
amplitude 0(6). The transverse spin (S+) t is the sum 
of an oscillating component cr° sc (<), a decaying compo- 
nent er^° c (i) with initial amplitude 0(6), and a secu- 
lar term a+°(t), which grows unbounded (linearly) in 
time. At fourth order in TL' V , (S z ) t also contains a sec- 
ular term. These difficulties, which have been reported 
previously^ii^ suggest the need for a more refined ap- 
proach. In the next subsection these problems will be 
resolved by working directly with the GME (in Born ap- 
proximation) to find the correct leading-order spin dy- 
namics for high magnetic fields. 
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B. Non-Markovian corrections 

Markovian dynamics are commonly assumed in spin 
systems j^L^I often leading to purely exponential relax- 
ation and decoherence times T\, and T%, respectively. For 
this reason, it is important to understand the nature of 
corrections to the standard Born-Markov approximation, 
and, as will be demonstrated in Section I VII on measure- 
ment, there are situations where the non-Markovian dy- 
namics arc dominant and observable. 



To apply the Born-Markov approximation to (S. 



+/t> 



change variables (S'l) t = e -H"»+")* (S+) t in (jUJ) and 
substitute E_)__|_(£) — > ESwt), which gives: 



(S'l) t = -iZ(S'i) t -i / dt'e-^'^Ut-niS'l)^ 

.In 

(52) 



We define the function ip(t) = 



where cu = u>„ + 

J t °° dt'e-^'^lit'), so that V(0) = ^++(s = iui). We 
find^a 



( S ») t = -i(m+z)(s';) t +ij t j q dt'iit-t')(si) t , 



The frequency shift ui is chosen to satisfy u> 



(53) 



-Re [V(0)] 



-Rc 



,(2) 



(s = z (w„ + uj)) 



to remove the 



oscillating part from (S") t . When \ui\ > 1, and after 
performing the continuum limit, we find a vanishing de- 



cay rate T 



Tm 



,(2) 



0. which shows that 



there is no decay in the Markovian solution for \lo\ > 1. 
After integrating the resulting equation, we have 



(s';) t = (s'i) +R + (t). 



(54) 



The Markovian solution is given by (S+) t = (S") , and 

the remainder term R+(t) = i J dt'ip(t — t') (S+) t , gives 
the exact correction to the Markovian dynamics (within 
the Born approximation). We rewrite the remainder term 



(55) 



R + (t) = i / dt'i)(t-t')((Sl) + R + (t')) 



Within the Born approximation, R+(t) is associated with 



a smallness 0(S 



N 



) (since ip{t) 



,(2) 



(t)), so the 



above expression can be iterated to evaluate the leading- 
order contribution to R+(t) in an asymptotic expansion 
for large w„. This gives 



R + (t)~-6I (c + + c-)(S + ) Q + e 



-iuj n t ^dec 
G , 



5 (t), (56) 



with af c {t) given in Eq. gTJ. 

Due to the inhomogeneous term N z (t) in l|29[l . the 
(S z ) t equation docs not have a simple convolution form, 
so it is not clear if a Markov approximation for (S z ) t 
is well-defined. However, applying the same procedure 



Symbol 



b' 



A 



Ho 



frif/Thf 



t c /r c 



Aq/2 = l,h 



b' + 2pIN 



N/uj„ 



ft 



f(c+ + c- 



Ap = A/N 



'ubB, — giu N B z 



b' + pi A 



A/2uJn 



2Nh/A 



N 
1-P 2 



S 2 = 



pIA 



l/(2p/) 2 iV 



2iVfi/yl 



Table I: Some symbols used in the text. The second column 
gives the value in dimensionless units, the third column gives 
the value in dimension-full units assuming An — ^, and the 
fourth column gives the value of each symbol in zero mag- 
netic field. The values shown are: the effective applied field 
b' , the total effective field (applied field and Overhauser field) 
seen by the electron cj„, the smallness parameter A, which 
determines the perturbative regime for electron spin dynam- 
ics, the smallness parameter S, which bounds the deviation of 
the electron spin from a Markovian solution, the coefficients 
c+ and c_ , in terms of the fraction of nuclear spins / = i up 
in the initial state /f , the electron spin precession frequency 
ila when the resonance condition uj„ = is satisfied, the time 
scale ihf for the decay of the electron spin in the presence of 
an initial /i z eigenstate of the nuclear system, and the time 
scale t c for the decay of the electron spin in the presence of 
an ensemble of initial nuclear spin states or a superposition 
of h z eigenstates at zeroth order in the nuclear spin-electron 
spin flip-flop terms. 



that was used on (5*+) t to determine the deviation of 
(S z } t from its initial value gives the remainder R z (t), to 
leading order in — , 



R z {t) ~ -2SI (c+ + c_) (S Z ) Q + 2pISI + af c {t). (57) 

Here, a z ec (t) is identical to the result from standard per- 
turbation theory, given by Eq. (|50|) . 

Corrections to the Markov approximation can indeed 
be bounded for all times to a negligible value by making 
the parameter 6 sufficiently small. However, the dynam- 
ics with amplitude O(S) are completely neglected within 
a Markov approximation. 

If wc use (S z ) t = (S z ) + R z (t) and Eq. JSU), and 
return to the rest frame for (S+) t , Eqs. I|56|) and l|57|) 
recover the high-field results from standard perturbation 
theory, given in Eqs. I|44|) and (|45|) . with one crucial 
difference. The result from standard perturbation theory 
contains a secular term, which is absent in the current 
case. Thus, by performing an expansion of the self-energy 
instead of the spin operators directly, the contributions 
that led to an unphysical divergence in (S+) t have been 
successfully re-summed. 





GaAs 


Si:P 


A 


90/ieF 


0.1 ^eV 


N 


10 5 


10 2 


B z 


7T 


0.1T 


P 








A 


0.25 


0.25 


S 


i(r 6 


10~ 3 




lOV 1 


10V 1 


Thf 


1 fiS 


1 j-LS 


T c 


5 ns 


100 ns 



Table II: Sample numerical values for the symbols listed in 
Table H] for a GaAs quantum dot or an electron trapped at a 
donor impurity in natural Si:P. 



C. Dependence on the wave function 

The purpose of this subsection is to evaluate the de- 
pendence of the non-Markovian dynamics on the form of 
the electron envelope wave function ip(r). The high- field 
dynamics, described by Eqs. 156fl and (|57f) . depend only 
on the integrals I± (i) . From Eq. H4U|) we find 



I±(t) 



dx I In x \ v xe ±ix \ 



d 
m 



1. 



(58) 



The time scale Thf for the initial decay of I± (t) is given by 
the inverse bandwidth (range of integration) of the above 
integral. In dimension-full units, Thf = -j*-. The long- 
time asymptotic behavior of I± (t) depends sensitively on 
the dimensionality d and the form of the envelope wave 
function through the ratio — . When — < 2, the major 
long time contribution to comes from the upper limit 
i«l corresponding to nuclear spins near the origin, and 
the asymptotic form of I± (t) shows slow oscillations with 

A ■ 



period 4lTh 



I±(t > 1) oc 



< 2. 



(59) 



When ■i- > 2, the major contribution comes from the 
lower limit x sa 0, i.e., nuclear spins far from the center, 
where the wave function is small. The resulting decay 
has a slowly-varying (non-oscillatory) envelope: 



I±(t > 1) oc 



ln u t 
t 2 ' 



1 > 1. 

to 



(60) 



Both of the above cases can be realized in physical sys- 
tems. For an electron with an s-type hydrogenic wave 
function bound, e.g., to a phosphorus donor impurity in 
Si, to = 1 and d = 3, which corresponds to the case in 
Eq. l|(jUfl . For an electron trapped in a parabolic quantum 
dot, the envelope wave function is a Gaussian (to = 2) 
and for d < 3, the asymptotics of I±(t) are described by 
Eq. (|59|) . These two cases are illustrated in Fig. [21 where 




t [2H/A ] 



Figure 2: Re [I+(t)/io] determined numerically from Eq. 1580 . 
For d — 3, m = 1 (solid line), this corresponds to a hydrogen- 
like s-type envelope wave function, and for d — m — 2 (dashed 
line), corresponding to a two-dimensional Gaussian envelope 
wave function. For the hydrogen-like wave function, nuclear 
spins far from the origin, with small coupling constants, are 
responsible for the slow (non-oscillatory) asymptotic behav- 
ior. In contrast, for the Gaussian envelope wave function 
nuclear spins near the center, with larger coupling constants, 
give rise to oscillations in the asymptotic behavior of I+(t)/Io. 



V. NON-MARKOVIAN DYNAMICS 

In this section we describe a complete calculation for 
the non-Markovian electron spin dynamics within the 
Born approximation. In the limit of a fully polarized 
initial state, our Born approximation applied to (S+) t 
recovers the exact solution of Ref. All results of this 
section are, however, valid for arbitrary polarization in 
high magnetic fields when the condition |A| <C 1 is sat- 
isfied. In addition, we find that the remainder term is 
bounded by the small parameter 6, \Rx(t)\ < 0(5), and 
the stationary limit (long-time average) of the spin can 
be determined with the much weaker condition 6 <C 1. 
In zero magnetic field, and for nuclear spin I = \, the 
relevant smallness parameter is S = (see Table QJ. 

We evaluate the Laplace transforms of (|29|l . 13Uf) : 
S x (s) = J °° dte~ st (Sx) t > Re[s] > 0, X = z, +, to con- 
vert the integro-differential equations into a pair of linear 
algebraic equations which can be solved to obtain 



S z (s) = 
S+(s) = 



(S z ) + N z (s) ^ 
s + iY, zz (s) 

(S+) 
s — iuj n + iS ++ (s) 



(61) 
(62) 



When the functions N z (s), Y, zz (s), S ++ (s) are known, 
the Laplace transforms in l|61|) and (|62|) can be inverted 
by evaluating the Bromwich contour integral: 



Re [/+(£) /io] is shown for d 



TO 



2 and d — 3, to = 1. 



7-M0O 



dse st S x (s), 



(63) 



7 — 200 
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the closed contour C shown in Fig. § c dse st Sx(s), 
gives 

(Sx) t + P x (t) = X = z,+, (64) 

i 

where the pole contribution Pi(t) = 
Res \e st Sx(s), s = s,-] is the residue from the pole 
at Si, and the branch cut contributions are 

0' (t) = - Im [e- iUnt Km] , (65) 

a=0,+,- 

= E ( 66 ) 

0=0, + ,- 

with branch cut integrals given by 



(67) 
(68) 



K z a {t) = / dse st S z (s-iu n ) 

K+{t) = f dse st S+(s). 
Jc„ 



Figure 3: The closed contour C used for evaluation of the 
inverse Laplace transforms of Sx{s), X = z, +. All non- 
analyticities of 1/D(s) are shown above, where D(s) is given 
in Eq. 17211 . Branch cuts are indicated by dashed lines, branch 
points by crosses, and open circles mark pole positions. The 
contour C a surrounds the branch cut extending from branch 
point 7 a = ai, a — 0,+,—. When the arc that closes the 
contour in the negative-real half-plane is extended to infinity, 
Cb becomes the Bromwich contour. The pole at S2 has finite 
real part and is present for b 7^ 0. The poles at si and S3 are 
always located on the imaginary axis. 



where all non-analyticitics of Sx(s) lie to the left of the 
line of integration. To simplify the calculation, here we 
specialize to the case of an electron confined to a two- 
dimensional parabolic quantum dot (d = m = 2), where 
the coupling constant integrals can be performed easily 
to obtain the explicit form for I±(s), given in Eq. (|42|l . 

Within the Born approximation, S z (s) has six branch 
points, located at ito n , i (u> n ± 1) , —iuj n , —i (uj n ± 1). 
We choose the principal branch for all logarithms, defined 
by log(z) = ln|z| + iarg(z), where — tt < arg(z) < n, in 
which case there are five poles in general. Three of these 
poles are located on the imaginary axis and two have fi- 
nite negative real part. S+(s) has three branch points (at 
s = 0, ±i), and three poles in general. One pole has finite 
negative real part and two are located on the imaginary 
axis. 

Applying the residue theorem to the integral around 



The contour C a runs from 7 Q — 00 + ir], around j a , and 
back to 7 Q — 00 — ir], where Tj — > + . The branch points 
are given by "f a = ai, a = 0, +, — , as illustrated in Fig. 
El In l|6l)|) we have used the fact that the branch cut in- 
tegrals for S z (s) come in complex conjugate pairs, since 
S z { s *) = [£z( s )]*- This relationship follows directly from 
the definition for the Laplace transform of the real quan- 
tity (S z ) t . 

Combining Eqs. Ip^. J33}> and (gSJ to obtain 

£zz (s — iu3 n ), and expanding in -3- gives 



(2) A 
iuj n ) = S ++ (s) + — (c + + C- 



0(6), (69) 



where we recall A = 



N 



and 6 = 



N 



The term 



— (c + + c_ ) gives rise to a small shift in the effective 
magnetic field experienced by (S z ) t . To simplify the pre- 
sentation, this shift is neglected, but it could easily be 
included by introducing a slight difference in the denom- 
inators of S z (s) and S+(s). This gives 

o / ■ , (S z ) + N?\ S -iLj n ) 
b z {s-iu n ) ~ y — , (70) 



D(s) 



S+(s) 



(S+) 
D(s) ■ 



(71) 



The denominator D(s) = s — it 
ator N?^ (s — iuj n ) are given explicitly by 



(2) 

iS\j_(s) and numer- 
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D(s) = s-ib' + Ns[c-log(s-i) + c + log(s + i)-(c+ + C-)log(s)}, (72) 
Ni 2 \s-iu n ) = -^( c+ + c _)-iA|[c + log(s + i)-c_log( S -i) + (c_-c + )log(s)]+0(<5). (73) 



The branch cuts and poles of S z (s — iu> n ) and S+(s), 
as given in Eqs. (|70fl and Q71[l. are shown in Fig. [3J We 
note that different analytic features will produce differ- 
ent types of dynamic behavior after the inversion integral 
has been evaluated. The branch cut contributions /3 X (t) 
have long-time tails that are non-exponential. Poles with 
finite negative real part will give rise to exponential de- 
cay. Poles on the imaginary axis away from the origin 
will lead to undamped oscillations, and a pole at the ori- 
gin will give a constant residue, independent of time. The 
rest of this section is divided accordingly, describing each 
type of contribution to the total time evolution of (Sx) t - 



A. Non-exponential decay 

The contribution to K x (t) circling each branch point 
7 Q is zero, so the branch cut integrals can be rewritten 
as 



K x (t)=e^ t dxe- xt i x {x, la ) 
Jo 



(74) 



where 



£x(x,j a )= lim \Sx(s%(x)+irj) 
i)— o+ 



Sx( S X (x)~i V )}, (75) 



with 



-X + 7 Q + 



—iu) n , X = z 
0, X = + 



(76) 



The form of K x (t) in Eq. I|74|) suggests a direct 
procedure for evaluating the long-time asymptotics of 
the branch cut contributions. For long times, the inte- 
grand of H74f> is cut off exponentially at x ~ j — ► 0. 
To find the asymptotic behavior, we find the leading 
x-dependence of £,x{x,"fa) for x — > + . We substitute 
this into (|74|) . and find the first term in an asymptotic 
expansion of the remaining integral. The leading-order 
long-time asymptotics obtained in this way for all branch 
cut integrals K x {t) are given explicitly in Appendix 
El When b' = 0, the denominator D{s) — ► when 
s — > 0, and the dominant asymptotic behavior comes 
from K x (t — > oo) oc A. For b' ^ 0, D(s) remains finite 
at the s = branch point and the dominant long-time 
contributions come from K x (t — > oo) oc - ^ - . In zero 
magnetic field, the leading-order term in the asymptotic 
expansion is dominant for times f > 1, but in a finite 



magnetic field, the leading term only dominates for times 
t > e\ b '\ /N . In summary, 



(3 x {t-> 1) oc 
x {f»e\ b '\l N ) oc 



1 

In? 
1 



ilnV 



b' = 0, 
b' ^ 0. 



(77) 
(78) 



This is in agreement with the exact result for a fully- 
polarized system of nuclear spins I = ^ in a two- 
dimensional quantum dot. This inverse logarithmic time 
dependence cannot be obtained from the high-field so- 
lutions of Section II VI The method used here to evalu- 
ate the asymptotics of the Born approximation therefore 
represents a nontrivial extension of the exact solution to 
a nuclear spin system of reduced polarization, but with 
|A| < 1 (see Table ID. 

The branch cut integrals can be evaluated for shorter 
times in a way that is asymptotically exact in a high mag- 
netic field. To do this, we expand the integrand of Eq. 
(|74|l to leading nontrivial order in — , taking care to ac- 
count for any singular contributions. For asymptotically 
large positive magnetic fields, we find (see Appendix lE|l : 



K x {t) ~ -i2nS {C x I+(t) + C x I-{t)) 



C 



x 



Net 



-zot 



(79) 



with coefficients C± given in l|51|l and in the above, 



za = x - ie(x ), 
x 



e(x) 



2irNc- 
x 

27rc„A ' Attc-X 



(80) 
(81) 

(82) 



In high magnetic fields, we will show that the exponential 
contribution to Eq. I|79|l cancels with the contribution 
from the pole at S2, P x {i)- We stress that this result 

is only true in the high-field limit ^J- ^> 1, where the 
asymptotics are valid. 



B. Exponential decay 

When b' = 0, there are no poles with finite real part. 
For b' ^ 0, a pole (at S2 in Fig. ^ emerges from the 
branch point at s = 0. The pole contribution P x (t) 



decays exponentially with rate 



-Re[s2], and has 
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-0.02 -0.01 0.01 0.02 
b'/N 



Although it does not correspond to the perturbative 
regime, it is interesting to consider the behavior of the ex- 
ponentially decaying pole contribution P 2 X (t) in the limit 
V — ► 0, since the Hamiltonian TL in Eq. is known to 
be integrable for B z — (b' = 0)m& For vanishing positive 
magnetic fields (b' — > + ), with logarithmic corrections in 
Jj-, where 6 = N(c+ + c_) and £ = exp {l + O (^)}: 



C&7&0 



W 2 



^2 + (0) 



(5+) 



iV( C+ + C _)ln(^)' 
{S Z ) Q - (c+ + c_) (2 pl 
N(c + + c.) In (& 



(89) 
(90) 
(91) 





b'/N 



where £ = c _^ c ■ The exponentially decaying contri- 
bution vanishes only when b' = 0, and does so in an 
interval that is logarithmically narrow. We have deter- 
mined the rate, frequency renormalization, and ampli- 
tude of the pole contribution P 2 + (i) numerically. The 
results are given in Fig. 0] along with the above asymp- 
totics for high magnetic fields, |6'| 3> N. 



Figure 4: Top: numerically determined rate F2 (solid line) 
and frequency renormalization U2 (dashed line) as a function 
of magnetic field b'/N. Bottom: iVRe [P 2 + (0)] (solid line) and 
iVIm [Pg (0)1 (dashed line) as a function of magnetic field for 
the initial state (S+} = (S x } = j. The dotted lines give 
the asymptotics for high magnetic fields from Eqs. 1851 , 18(i> . 
and JH3- The parameters used were p = 0.6, N — 10 , 1 = 1 



an envelope that oscillates at a frequency determined by 
lj 2 = Im[s 2 ]: 



2 



e -r a t e -t(a,„-u, a )t p *( ) 



-V2t lLJ2t p-f - 



^2 + (0). 



(83) 
(84) 



Setting S2 = —T2 + i^2i wc nn d the decay rate T2, fre- 
quency renormalization L02 , and amplitudes of these pole 
contributions from asymptotic solutions to the pair of 
equations Re[D(s 2 )] = lm[D(s 2 )} = and P 2 A "(0) = 
Rcs[Sx(s),s = S2] for high and low magnetic fields b' . 
T 2 , u>2, and {0) have the asymptotic field dependences 
(for high magnetic fields b' N): 



To 



CJ 2 



±- 



" 2?r Nc^ 



ui n ^ 0, 



±- 



± 



P2 X (0) 



~2ttc t N 
1 =F i2nNc T 



47rc T r 2 



£ 0, 



LU n ^ 0. 



(85) 
(86) 
(87) 



C. Undamped oscillations 

The point si in Fig. [3] corresponds to s — for S z (s), 
so undamped oscillations in (S z ) t arise only from the pole 
at s 3 : 



Pg(t) = e- l(w ™- W3)t P 3 z (0). 



(92) 



Both poles on the imaginary axis give undamped oscilla- 
tions in 

P+{t) + P+(t) = e IWlt P+(0) + e^ 3t P 3 + (0). (93) 
For high magnetic fields, \b' /N\ ^> 1, 

w 1/3 ~ b' + 2pIN = u n) b' ^ 0, (94) 

W3/1 ~ T-l T f± exp (-^) , (95) 

where /± = (§)(^±) (l + 0(^)). The frequency in 
Eq. (|94|l corresponds to a simple precession of the elec- 
tron spin in the sum of the magnetic and Overhauser 
fields. The second frequency, Eq. I|95(l . describes the 
back-action of the electron spin, in response to the slow 
precession of the nuclear spins in the effective field of the 
electron. 

For large the pole corresponding to simple preces- 
sion is dominant, while the other has a residue that van- 
ishes exponentially: 
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10 20 30 40 



t [2N/T/A] 

Figure 5: Contributions to the inverse Laplace transform of 
{S z ) t . We show the envelopes of the rapidly oscillating func- 
tions 2NRe [P£(t)] + 0.3, 2iVRe [Pi{i)\ + 0.2, and -Nf3 z {t), 
determined numerically. The long-time asymptotics of f3 z (t) 
from Appendix El are also shown (dashed line). The sum 
of all contributions is used to obtain the population of the 
spin-up state: P](t) — \ + {S z ) t (inset). The electron spin 
begins down: {S~) Q = —\. Other parameters were I = |, 
N = 10 5 , b' = f (this value of b' gives, e.g., B z ~ IT in 
GaAs), and p = 0.6. The time t is given in units of = ^S. 
for d = m = 2 in Eq. fgJJ (2^ ~ lps in GaAs). These 
values correspond to the weakly perturbative regime, with 
A = Yj < 1. Note that p T (t) < i for all times. 



Within the Born approximation, we find 



(S z ) +p/<S + O 






l + (c + +c_)5 + 


(ft) 



(100) 



The result in Eq. (|TUU)l follows from Eqs. (JHTJ, J2TJ|, 
(I32[) . (|34|) . JsHJl, and (|37|> by expanding the numera- 
tor and denominator in using the coupling con- 
stants Afc = 2e~ k ' N and performing the continuum limit. 
(5 Z ) gives the stationary level populations for spin-up 
and spin-down: py^ = \ ± (5 l z ) oc , which would be fixed 
by the initial conditions in the absence of the hypcrfine 
interaction. This difference in p^ ^ from the initial val- 
ues can be regarded as leakage due to the nuclear spin 
environment. We note that the stationary value depends 
on the initial value (S z ) , from which it deviates only by 
a small amount of order 5. This means, in particular, 
that the syst em is non-ergodic. We will find that correc- 
tions to (>S' Z ) 00 at fourth order in the flip-flop terms will 
be of order S 2 , so that the stationary limit can be deter- 
mined even outside of the perturbative regime |A| < 1, 
in zero magnetic field, where 6 = for I = h, provided 

P»VN- 

E. Summary 



A + /3 (o) 

K/i (o) 
Pg(0) 



l + ±(c++c_)<5' 



(96) 



(S. 



2 

+; 



Nc± 
b 1 



/±exp 



\b'\ 



c±N 



, b'^0, (97) 



^exp^-J^L), y >0 . (98 ) 



When the magnetic field b' compensates the nuclear 
Overhauser field [^ 2 ] n „ (w n ~ 0, the usual ESR resonance 
condition in the rotating frame), the poles at points S\ 
and S3 have equal weight, and arc the dominant contribu- 
tion to the electron spin dynamics. Since the resonance 
condition corresponds to the strongly non-perturbativc 
regime, |A| S> 1, we delay a detailed discussion of the 
resonance until Section IVlIJ 



D. Stationary limit 



The contribution to (S z ) t from the pole at s — gives 
the long-time average value (S z ) , which we define as 
the stationary limit: 



(S*)oo = l™ i / (S,) t dt = lim sS z (s). (99) 



The results of this section for low magnetic fields are 
summarized in Fig. [5] which corresponds to the weakly 
perturbative case, |A| < 1, and displays all of the dy- 
namical features outlined here. 

In very high magnetic fields (b' N), corresponding 
to the strongly perturbative case, we combine Eqs. I|79(l . 
(|MJl . 0, and l|100|) to obtain the asymptotic forms to 
leading order in — : 



(S+) t 
{S,) t 



: (*) 



T dcc 



dec 



(*), 
(*), 



(101) 
(102) 



where the functions cr^ sc (0, crf c {t), (S z ) 



_dcc 



(*), 



given in Eqs. |@BJ|, l|17 )l .l|39 ^l . and (JSUJl are evaluated 
for d = m = 2. We stress that <Jx c (t) oc S -C 1 is a 
small fraction of the total spin. The exponentially decay- 
ing contribution from P% (t) is canceled by the exponen- 
tial part of the high- field branch cut, given in Eq. J75J). 
This result is in agreement with the high-field asymptotic 
forms found earlier in Section llVI Numerical results for 
the level populations pj 1^ (t) = ^ ± (S z ) t are given in Fig. 
[H] along with the above asymptotic forms. The secular 
term that appeared at lowest order in the standard per- 
turbation expansion of (S+) t is again absent from the re- 
sult obtained here via the GME. At fourth order, f-linear 
terms also appear in the standard perturbation expansion 
for the longitudinal spin (S z ) t M^ Due to the numera- 
tor term N z (s) in the expression for S z (s) (Eq. i|61|) ). 
it is not clear if all divergences have been re-summed 
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t [2NH/A = 1 us] 

Figure 6: Envelope of the time-dependent spin level popu- 
lations in high magnetic fields. We give results from nu- 
merical inversion of the Laplace transform (solid line) and 
the asymptotic branch cut integral for high magnetic fields 
combined with numerical results for the pole positions and 
residues (dashed line) . Top: spin-down level population when 
the electron begins in the up state, along the nuclear spin po- 
larization direction ({S z ) = §)■ Bottom: spin-up population 
for an electron that begins pointing in the opposite direction 
((S z ) = — |)- The parameters used were N = 10 5 , p = 0.6, 
1=2 an d b' — 8JV, corresponding to a field of B z = 14 T in 
GaAs. 



an initial ^-pulsc. Each spin precesses in its own local 
effective magnetic field w„. The phase factor e 1 "™* winds 
in the "forward" direction for a time r. The sign of ui n 
(direction of the local magnetic field) is then effectively 
reversed with a 7r-pulse along the x-axis: u n — > — uj n . 
The phase factor e~ lUnt unwinds in the following time 
interval r, and the electron spin magnetization rcfocuses 
to give an echo when the phase factor e ~ % ^ n2T = 1 si- 
multaneously for all spins in the ensemble. As is usually 
assumed, we take the pulse times and measurement time 
during the echo to be negligible^ The spin echo enve- 
lope gives the ensemble magnetization (the electron spin 
expectation value) at the time of the echo as a function of 
the free evolution time 2r before the echo. We note that 
the decaying fraction of (S z ) t , crf ec (i), also precesses with 
the phase factor e tuJnt (see Eq. ||SDJ|), so the same pulse 
sequence can also be applied to measure the decay of the 
longitudinal spin, omitting the initial -|-pulse. The Hahn 
echo envelope should show a small initial decay by 0(S) in 
a time scale Thf due to the contact hyperfine interaction, 
followed by a slow decay due to spectral diffusionSi*^ 2 ^ 
with a time scale Tdd ~ 1CP 4 s. We note that a rapid ini- 
tial decay of the Hahn echo envelope has been measured 
for natural Si : P, but is absent in isotopically enriched 
28 Si : P, in which no nuclei carry spin4^ 



for (S z ) t in the pcrturbative expansion of the self-energy. 
This question is addressed in Sec. IVIII with an explicit 
calculation of the fourth-order spin dynamics. 

In the next section we propose a method that could be 
used to probe the non-Markovian electron spin dynamics 
cxp er imcnt ally. 



VI. MEASUREMENT 

In high magnetic fields (&' S> N), the decaying fraction 
of the electron spin is very small (O (6 rts 57?)) ■ Never- 
theless, the large separation between the hyperfine in- 
teraction decay time (rhf = ~ l/^s) and the dipolar 
correlation time (Tdd ~ 100 /is in GaAs) of the nuclear 
spins should allow one to obtain valuable information 
about the electron spin decay from a conventional spin 
echo technique applied to an ensemble of electron spins. 

In principle, the non-Markovian electron spin dynam- 
ics should be visible in the electron spin echo enve- 
lope obtained by applying the conventional Hahn echo 
sequenced ^ — t — ir x — r — ECHO to a large ensemble of 
electron spins. This can be done by conventional means 
for an electron trapped at donor impurities in a solidi^ 
or from a measurement of transport current through a 
quantum dot^i^ The effect of this echo sequence can 
be summarized as follows. The electron spins are ini- 
tially aligned along the external magnetic field B z . At 
time t = the spins are tipped into the x — y plane with 



The fraction of the spin that decays in the time Thf 
is small, of order 5, in the perturbative regime. It 
may be difficult to detect this small fraction using the 
conventional Hahn echo. This problem can be reduced 
by taking advantage of the quantum Zeno effect, using 
the Carr-Purcell-Meiboom-Gill (CPMG) echo sequence 
§ - (r - tt x - r - ECHO - r - tt^ - t - ECHO) rcpeat . 
During each free evolution time between echoes, the elec- 
tron decays by an amount of order d. At each echo, a 
measurement of the electron spin magnetization is per- 
formed. For a large ensemble of electron spins, this mea- 
surement determines the state ps of the electron spin 
ensemble, forcing the total system into a direct product 
of electron and nuclear states, as in Eq. (US) . Repeti- 
tion of such measurement cycles will then reveal the spin 
decay due to the hyperfine interaction (by order 6 af- 
ter each measurement) until the magnetization envelope 
reaches its stationary value. If the electron spin decays 
during the free evolution time due to spectral diffusion 
with a Gaussian envelope, then we require the condi- 
tion (^r^j "C 8 <C 1 for the effect of spectral diffu- 
sion to be negligible compared to the effect of the hyper- 
fine interaction^ The non-Markovian remainder term 
gives the total change in electron spin that has occurred 
during the free evolution time 2r: Rx (2T)| e ±i„„2-r_i = 
{S x ) 2T -{S x ) a \ e±Wn2T=1 = M x (2t) - M x (0), where 
Mx(t) is the CPMG magnetization envelope. In high 
magnetic fields, and when there are many echoes before 
the magnetization envelope decays, the CPMG magneti- 
zation envelopes M x (t) will therefore obey the differcn- 
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Figure 7: Longitudinal decay rate — jy of the CPMG echo 

J i 

envelope as a function of the free evolution time 2r between n- 
pulses for an electron trapped at a phosphorus donor impurity 
in Si:P (top) and in a two-dimensional GaAs quantum dot 
(bottom). The free evolution time is given in units of « 

(the equality is exact for d = m = 2 in Eq. Ij43^ ) . In 
a GaAs quantum dot containing N = 10 5 nuclei or for an 
electron trapped at a shallow donor impurity in Si:P with 
N — 100 nuclear spins within one Bohr radius, ~ \ 
We have used I = |, p = 0.6, and magnetic field values from 
Table ITTlto determine the frequency units on the vertical axis. 



tial equations 



dt 



M x {t) 



Rx(2t) 



2t 



X 



(S x ) =M x (t),ei 



(103) 

where the high- field expressions for Rx(t), given in Eqs. 
(|56|l and i|57|) . should be used. Thus, the decay rate of 
the CPMG echo envelope Mx, as a function of the free 
evolution time 2t, is a direct probe of the non-Mar kovian 
remainder term Rx(t). 

Since the magnetization envelopes Mx(t) are found as 
the result of an ensemble measurement, it is necessary to 
perform an average over different nuclear initial states \n) 
that may enter into the solutions to Eq. (|103fl . The local 
field-dependent phase factors have been removed by the 
echo sequence, so the only effect of the ensemble average 
is to average over S = and c±, which appear in the 

overall amplitude of -^Mx(t). The relative fluctuations 
in these quantities are always suppressed by the factor 



'N 



for a large nuclear spin system. 



In the high- field limit, we find the longitudinal and 
transverse magnetization envelopes M z (t) and M + (t) de- 
cay exponentially with time constants T± and T^ 1 = 
2T-j M , respectively. M+(t) decays to zero, and M z (t) de- 
cays to the limiting value 



M z (oo) 



lc_ 



pi 



2 c_ + c + c_ + c + 
For nuclear spin I = 4, M z (po) = i.e., the electron 



magnetization acquires the polarization of the nuclear 

spin bath. However, since c± cx I 2 , M z (oo) — > in the 

large-spin limit. Thus, a larger fraction of the electron 

spin decays in the limit of large nuclear spin. We give 

plots of the longitudinal spin decay rate for M z (t), tJw, 

1 1 

as a function of the free evolution time 2r for two types of 
envelope wave function in Fig. [7| These plots have been 
determined by integrating Eq. (|103|l using the high-field 
expression for R z (t) given in Eq. 1)57(1. No ensemble aver- 
aging has been performed to generate these plots. When 
2t <C Thf, the envelope decay rate increases as a function 
of 2r as more of the electron spin is allowed to decay be- 
fore each measurement. The rates reach a maximum at 
some time 2r w Thf, and for 2t 3> Thf, the electron spin 
saturates at its stationary value and the envelope decay 
rates oc i are determined only by the free evolution 
time. Note that there are slow oscillations in the CPMG 
decay rate for an electron in a GaAs quantum dot, with a 
Gaussian wave function, but none for an electron trapped 
at a donor impurity in Si: P. 



VII. BEYOND BORN 

The goal of this section is to address the range of va- 
lidity of the results obtained in Sec. First, we show 
that the Born approximation for (S+) t recovers the exact 
solution for I = ^, p = 1. We then discuss the behav- 
ior of the Born approximation near the ESR resonance, 
where uj n ps 0. Finally, we consider the expression for 
(S z ) t , obtained by including all fourth-order corrections 
to the reduced self-energy, and show that our expression 
is wcll-bchavcd in the continuum limit. 



A. Recovery of the exact solution 



When I = | and p = 1, we have c_ 

(2) 



1 and c + = 0, 



which gives SV+( S ) = -JE 



from Eq. 1A22 



We insert this into l|rj2"|) and use w n = b' + | J2k ^-k = 



V + 4 to obtain 



S+(s) 



(104) 



A k /2 



The Schrodingcr equation for a state of the 
format)) = Ofr(t)|ftTT-"> + ^(t)|W-) + 
J2k Pk{t) |1tT ■ ' ' ik] ■ ■ ■), where the large arrow gives 
the state of the electron spin and the thin arrows give 
the states of the nuclear spins, has been written and 
solved (for a fully polarized nuclear spin initial state, 
j3k(t = 0) = V k) in Laplace space to find the long-time 
asymptotic electron spin dynamics previously^ In Ref. 
Ha the symbol a(t) was used in place of aqj.(f). The 
fully-polarized state Iff "ft • • ■ ) is an eigenstate of the full 
Hamiltonian T~C , so a^(0), which 



allows us to write S+(s) = a^(t = 0)a^ (s — § (b 1 + =§)). 
We solve the time-dependent Schrodinger equation for 
\tp(t)) in Laplace space, giving 



"4 ( s ') 



a^(t = 0+) 



i(y + 4) + jE 



(105) 



where s' = s — | (6' + Thus, in the limit of full po- 
larization of the nuclear system, the Born approximation 
applied to (S+) t becomes exact. For a fully polarized 
nuclear spin system (S z ) t is given by the relationship 

2^ 



l(l-2Mi)l) 



1 



(S,. 



4(t=o) 



Un- 



fortunately, this result is not recovered directly from the 
Born approximation for (S z ) t , as we will show in the next 
subsection. 



B. Resonance 

On resonance, uj n = 0, i.e., the external field 6' com- 
pensates the Ovcrhauscr field [/i 2 ] nn . The resonance 
is well outside of the pcrturbative regime, defined by 



IA 



N 



< 1, but we proceed in the hope that the Born 

approximation applied to the self-energy captures some 
of the correct behavior in the non-perturbative limit. On 
resonance, the major contributions to (S z ) t come from 
three poles, at s = 0, s = S3, and s = s^: 



(Sz) t *(S z ) 00 + 2Re[Pi(t)] 



(106) 



Before applying the continuum limit, the stationary limit 
for (S z ) t is 



(S, 



(S z ) + \(c.-c + )N tot 



1 



N 



(107) 



2 - 1 ' tot 

After applying the continuum limit, 7V to t — ► 00, we obtain 
1 c_ - c+ pi 



2 c. 



c_ + c + 



(108) 



i, (5' 2 ) 00 = |, which appears to be an intuitive 



For I 

result. However, evaluating the remaining pole contri 
butions at the resonance, we find, for a two-dimensional 
quantum dot, 



2Re[P 3 (t)} = 
where 



(S z 



2pl 



C- + C 1 



cos (Sl t) + O 



N 



+ c_). 



1 

N 
(109) 



(110) 



The results in i|108|) and (|109|l do not reproduce the exact 
solution in the limit p = 1, I = i, and do not recover the 
correct t = value of (S z ) t . The Born approximation 




B -B" [T] 



Figure 8: {S z ) x evaluated within Born approximation near 
the resonance, from Eq. ill ( lilt where = — 2g ^ B ■ We 
have used the value of j4 for GaAs, g* = —0.44, = 10 5 , 
and / = i. (S l ») = — ^ for all three curves and results 



are given for p — (solid line), p — A (dotted line), and 
p = 1 (dashed line). The vertical dash-dotted lines indicate 
the magnetic fields where the relevant smallness parameter is 
unity: \S\ = 1. 



for (S z ) t , as it has been defined here, breaks down in the 
strongly non-perturbative limit, although the transverse 
components are better behaved. 

On resonance, the poles at s\ and S3 are equidistant 
from the origin, and the major contributions to (S+) t 
come from these two poles: (S+) t w Pi(t) + P^t). Eval- 
uating the residues at these poles, 



(S. 



+/t 



(S 



l-O 



cos(n t), (111) 



which suggests that a fraction O (jj) of the spin un- 
dergoes decay, and the rest precesses at a frequency 



Do- When / 



and in proper energy units we have 



/SN 



from Eq. (|110jl . While it does not violate 



positivity, as in the case of (S z ) t , this expression should 
not be taken seriously in general, since this result has 
been obtained well outside of the pcrturbative regime. 
The above does, however, recover the exact solution in 
the limit p = 1. We show the stationary limit of (S z ) t in 
Fig. |H1 using typical values for an electron confined to a 
GaAs quantum dot. 



C. Fourth-order corrections 

The fourth order expansion of the self-energy for (S z ) t 
is given in Appendix lAl The discrete expression for the 
numerator term Ni 4 \s) contains second order poles (sec- 
ular terms). The fourth-order expression for S z (s) inher- 
its these second order poles (see Eq. (IfTlT ). When the 
Laplace transform is inverted, this will result in pole con- 
tributions that grow linearly in time. However, when the 
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continuum limit is performed, which is strictly valid for 
times shorter than t « \Tn (see Appendix EJ , all poles 

in Ni 4 \s) are replaced by branch cuts. The integrals 
around the branch cuts can then be performed to obtain 
a solution for (S z ) t , valid for times t < s/N. 

All relevant non-analytic features (branch points and 
poles) of S z (s) occur in two regions of the complex plane: 
about the origin s 0, and at high frequencies, around 



s w ±iuj n . Inserting an initial nuclear state \n) for a large 
uniform system (see Appendix IB")) , expanding the fourth- 
order self-energy to leading order in -j- about the points 
s = and s = —iuj„, performing the continuum limit, 
and evaluating the integrals over coupling constants, we 
obtain (where the overbar and "conj." indicate complex 
conjugate for s real): 



A^ 4) (s - iu> n ) ~ -— {c+c_ [Li(a) + L 2 {s) - L 3 (s) - conj.] + ^(s) - c^L^s)} 
EW(* " ^n) - -NA {c+c_ [Lj (a) + L 2 {s) - L 3 (s) + conj.] + c^L^s) + cLL^a)} 



isS 2 



3(piy 



14c + c_) s 2 L A {s) 



with coupling constant integrals Lj(s) given by 



Li(s) 



2(s + i) 2 



[log(s + i) - log(s)] 



L 3 (s) = [slog(s + z)-slog(s)-z] 2 

L 4 (s) = i - i [s 3 + 3s + 2i] pog(s + i) - log(s)] - i [s 3 + 3s - 2i] [log(s - i) - log(s)] 
6 6s 6s 



(112) 
(113) 
(114) 

(115) 



(116) 

(117) 
(118) 



and 



L 2 (s) = log(s + i) — log(s) - i [(s + i) log(s + i) - (s + 2i) log(s + 2i) + s log(s) - (s - i) log(s - i)] 



s+'i 



du 



log(2u — s — i) — log(2u — s) 



(119) 



Noting that lim s ^o s 2 ^4(s) = 0, we find the corrections to the stationary limit for (S z ) t . At fourth order in the 
flip-flop terms, this gives 



{ S z ) 0+P I5+l(cl-cl)5 2 + o(%) 
1 + (c+ + c_)<5 - 3( P I) 2 S 2 + O fe) 
I 



(120) 



The fourth-order corrections to the self-energy at high 
frequency (s m —iLu n ) are suppressed relative to the Born 
approximation by an additional factor of the smallness 
parameter A, as expected from the analysis given in Ap- 
pendix ^ However, the low-frequency (s w 0) part of 
the fourth-order self-energy is suppressed by the much 
smaller parameter 5. This allows us to determine the 
stationary limit of (S z ) t with confidence even when the 
magnetic field is small or zero, provided the polariza- 
tion is sufficiently large. When b' = and / = |, we 
have S = -^jj , so the stationary limit can be determined 
whenever p 3> — j= . 



It is relatively straightforward to find the time- 
dependence as t — > oo for the S z branch cut integrals at 
fourth order. Neglecting contributions from the branch 
cuts near s ~ 0, which arc suppressed by the factor 
S 2 , and when p < 1 so that the coefficient c + ^ 
(c.f. Eq. (J51J), we find the major contributions at long 
times come from the branch points at s = ±i, where 
L 2 (s) oc log (s + i). For any magnetic field, we find: 



R z (t — > oo) oc 



1 



t In' 3 1 



(121) 



For 6' ^> 2pIN, this time-dependence will be domi- 
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L S ~ E s 


E s ~ E<?> 


E s ~ E« 


Es^E^+E^ 


Es * Ef 


E S ~ E<?> 




b' = 0, 


b'^0, 


6' t^O, 




|A| « 1, 


|A| « 1, 




- = 1, 


- = 1, 


- = 1, 


£ = 1, 


^<2, 


± > 2. 

m — ' 




l 


PJ 1 £ i > 1 


t » e M /Ar » rj 1 


£ » el b l /JV , 6' > 2p/iV 


t > 1 


1 


Rxit) oc 


l/lnt 


e ^ 2 t e -r 2 t 


l/iln 2 i 


1/Hn 3 = 2 


(l/i)£ e ±" 


ln^f/* 2 , v=--l 

1 ' m 



Table III: Results for the decaying fraction of the spin (|7?x(i)| < 0(5) Vi) in various parameter regimes. Results are given for 
both remainder terms Rxit), X = z,+, within the Born approximation for the self-energy Eg ~ Eg 2 ' and for R z (i) at fourth 
order in the nuclear spin-electron spin flip-flop terms Es ~ Eg 2 ' + Eg 1 ' when p / 1. The first three columns are exact in the 
limit of full polarization (p = 1) of the nuclear spin system, but still may describe the correct electron spin dynamics in the 
weakly perturbative regime, |A| < 1. The last two columns give the correct electron spin dynamics in the strongly perturbative 
regime, |Aj < 1. 



nant when t exp 



Thus, wc find that the 



§(£) J) to ^ 



fourth-order result has a faster long-time decay than the 
Born approximation, and that the associated asymptotics 
are valid at the same times as the Born approximation 
asymptotics (see Eq. <|78|) ). Thus, higher-order correc- 
tions may change the character of the long-time decay in 
the weakly perturbative regime, where they are not neg- 
ligible. In contrast, in the strongly perturbative regime 
|A| <C 1, the fourth- and higher-order terms are negligi- 
ble, so the Born approximation dominates for all times 
t < exp(\b'\/N). 



VIII. CONCLUSIONS 

We have given a complete analytical description for the 
dynamics of an electron spin interacting with a nuclear 
spin environment via the Fermi contact hyperfine inter- 
action. In a large magnetic field, our calculation applies 
to a nuclear spin system of arbitrary polarization p and 
arbitrary spin /, prepared in an eigenstate of the total 
^-component of the (quantum) nuclear Ovcrhauser field. 
In the limit of full polarization p = 1 and nuclear spin 
I = \, the Born approximation applied to the self-energy 
recovers the exact dynamics for (S+) t and (S z ) t , with all 
non-perturbative effects. We have shown explicitly that 
the dynamical behavior we calculate in Born approxi- 
mation is purely non-Mar kovian, and can be obtained 
in the limit of high magnetic fields directly from the re- 
mainder term to a Born-Markov approximation. By per- 
forming our expansion on the self-energy superoperator, 
we have re-summed secular divergences that are present 
in standard perturbation theory at lowest (second) or- 
der for the transverse components (S+) t and at fourth 
and higher order for the longitudinal spin (S z ) t . For low 
magnetic fields b' < N , but still within the perturbative 
regime (|A| < 1), the Born approximation for the elec- 
tron spin shows rich dynamics including non-exponential 
(inverse logarithm) decay, exponential decay, and un- 
damped oscillations. For high magnetic fields V 3> N, 
and for — < 2, the electron spin shows a power-law de- 
— 

cay (~ (t) m m d-dimensions for an isotropic envelope 



wave function of the form ip{r) oc exp 

stationary value with a time scale ~ ~Spi m agree- 
ment with the exact solution for a fully polarized nuclear 
spin system^i*i£ Above a critical ratio, ^ > 2, the spin 
decay asymptotics undergo an abrupt change, signaled 
by a disappearance of slow oscillations in the decay enve- 
lope. We have summarized these results in Table ITTT1 Wc 
have also suggested a method that could be used to probe 
the non-Markovian electron spin dynamics directly, us- 
ing a standard spin-echo technique. We emphasize that 
the electron spin only decays by some small fraction of 
its initial value, of order 5 (see Tables Q] [nj), an d the 
decay is generically non-exponential at long times (see 
Table 1111(1 . The results of this work may therefore be of 
central importance to the development of future quan- 
tum error correction schemes, which typically assume an 
exponential decay to zero. The fact that the stationary 
value of the spin depends on the initial value implies that 
this system is non-ergodic. Based on this observation, we 
postulate a general principle, that non-ergodic quantum 
systems can preserve phase-coherence to a higher degree 
than systems with ergodic behavior. It would be inter- 
esting to explore this connection further. 
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Appendix A: SELF-ENERGY EXPANSION 

To expand the self-energy superoperator E5 in powers 
of Ly, we have found it convenient to work in terms of 
a superoperator matrix representation. Here wc give a 
brief description of its use and apply it to generate the 
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reduced self-energy at second order in Ly for all spin 
components, and the fourth order for the longitudinal 
spin. 

Any operator O that acts on both the electron spin and 
nuclear spin Hilbert spaces can be written in terms of the 
2x2 identity, oo, and the Pauli matrices erj, i = [x, y, z): 



O = ^2 CiCJi 

i— (0,a:,y ,z) 



(Al) 



where the coefficients c, are operators that act only on 
the nuclear spin space. Equivalcntly, O can be writ- 
ten in terms of the operators p^n = | (cr ±ff 2 ), S± = 
\ (a x ± ioy), i.e.: 



O = fc T p T + k i p i + fc+SL + k-S + 



(A2) 



with operators kj that act on the nuclear spin space. We 
have labeled the coefficients kj in this way so that when 
O = ps is the electron spin density operator, k± = (S±). 
A superoperator S acting on O maps it to the operator 
O' with new coefficients: 

SO = O' = k\ P] + k[ Pl + k' + S- + k'_S+. (A3) 

This allows us to write O as a vector and iS as a 4 x 4 
matrix, the elements of which are superoperators that act 
on the nuclear spin space, and are determined by (|A2() 
and Eat: 



O = (k^k^k^k.) 1 



(A4) 
(A5) 



C = [5] and fc^ = ^ ScK/sfys, where a, /3 =|, j, +, -. 

Laplace transforming the reduced self-energy given in 
(HJ yields 



^s(s) = -iTrji- 



1 



s + 

which is expanded in powers of Ly 



L VPI (0), (A6) 



1 



iQL 



fc=0 



1 



iQL Q 



—iQL\ 



1 



s + iQLc 



(A7) 



To obtain these higher order terms in the self-energy, 
we form products of the free propagator s+i g La and the 
perturbation QLy. The free propagator is diagonal in 
the basis of {pyi, S±}, and is given in terms of 2 x 2 
blocks by 



iQL 




Gf(s) 



where 




G°(s) 

G°(s) 

G° + (s) 

G°_(s) 



(A8) 

(A9) 
(A10) 



In the above, 



1 



G°±(s) 



1 



(All) 
(A12) 



where we define the new (nuclear spin) Liouvillians 
by their action on an arbitrary operator O: L^O = 
[cj, 0] ± , lo = b' + h z . The perturbation term contains 
only off-diagonal elements when written in terms of 2 x 2 
blocks: 



[QLy] = 



Vj 
V' I 



(A13) 



where we find 



Vr = ^ 



2 V -hi 



Q h% -h<> 
1 2 1 -h R h L _ 



(A14) 

In the above expression, we have introduced superoper- 
ators for right and left multiplication: 



O r A 
O h A 



AO 
OA. 



(A15) 
(A16) 



Only even powers of Ly can contribute to the final 
trace over the nuclear system, so we consider a general 
term in the expansion of the self-energy 



[QL X 



1 



iQL 



2 k- 



Sfc o 

o £1 



Z k = (VjGfVlG ^ 



4 = (VlGVjViGff , 



(A17) 



(A18) 



By inspection of the form of V/, Vj, we find that the 2x2 
matrix TrjY,' k pj(0) is diagonal when p/(0) = \n) (n\, and 
\n) is an eigenstate of h z (as in Eq. I|20(l ). since the off- 
diagonal components always contain terms proportional 
to h\ or h 2 _. Thus, to all orders in the perturbation Ly, 
the reduced self-energy takes the form 



( S TT ( S ) £ n (s) 
E iT («) S u (s) 

£ ++ (s) 




V 














\ 



(A19) 



-(«)/ 



The number of matrix elements left to calculate can be 
further reduced with the relationships S|| (s) = — Sj^ (s), 
^Tl( s ) = — ^li ( s ) i which follow directly from the con- 
dition Trps = => p T (i) = -f>i(t) and the GME 

Pa = -*'S/3=T4 jl dt '^ap{t - f)pp(t'), a =|,|- % di- 
rect calculation we find 
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sg>(.) 

4 2 i(*) 



i 



S - lLU n k S + ILUnk 
1 1 



\ ^2(i h +}nk i h -}kn + \h-\nk \ h +] kn ) ~ _ 



In the above, uj n k = \ (w n + a;*,), buj nk = § (w„ — Wfe), and a/j = 6' + [/izLj- At fourth order, 



(A20) 
(A21) 
(A22) 



^1"T ( s ) 



16 



E ^"U [h + ) k3k2 [h-] k2kl [h + ] kin (1 - *„ fe2 K^ fe 3 (s) + a 4 fe ^( S 



. k\k 2 k 3 



\ E [ft-U [*+]*,« i h +U \h-\ kin > (A23) 



4f( s ) = { E t^+U [*»-]*.*, [fc-iM ((i - <w^ fe2fe3 (*) + °lt k3 {s)) 

~ E [^U [ ft -lfaT. [*-W [^ + ] fcin ^c fc2 ( 5 )} , (A24) 

k±k 2 ) 



where the overbar indicates complex conjugation for s real and 

fclfc 2 fc3 ' ~ ' 



s - i5u nk2 s - iu nkl \s — iui nk3 s + iuj k2k3 



kxk 2 k 3 
'4B 



T kik 2 
'AC 



(s) = 



» 



1 



1 



s - iSu} kl k 3 \s - iu klk2 



Every two powers of the perturbation Ly are associ- 
ated with an additional sum over sw N nuclear spin sites, 
since every spin flip up must be paired with a flop down. 
Non-analyticities (poles) of the self-energy occur in two 
regions of the complex plane: at high frequencies, near 
s f=a ±iuj n , and at low frequency, around s « 0. Expand- 
ing near either of these two points gives an extra factor 
£P for every two orders of QLy s _ 1 1 q Lq ■ The self-energy 

at (2k) th order is then suppressed at least by the factor 
A fc . where A = — : 



E^(s) cx A A 



E^ 2fe) (s - iw n ) oc A* 



(A28) 
(A29) 



1 



1 



1 



1 



s + iSu nk2 s + iuj nk3 \s + iuj nkl s - iu; klk2 



1 



1 



1 



S + lU) k2 k 3 J \S + lU) n h 3 s - lLO kl71 



r 



, (A25) 

(A26) 
(A27) 



Thus, in general, for the perturbation series to be well- 
controlled, we require |A| <C 1. 



Appendix B: COEFFICIENTS c± 

We are interested in evaluating the expressions given in 
Eqs. HA20|> . (|A21|) . and (|A22jl . To do this, we investigate 
objects of the form 



E^U [M** /*(*). 

k 



(Bl) 
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where fe(k) is a function of the state index k. Inserting 
| re), as given in Eq. (|20p. into (|B1|I . we find 

E [ ft ±U IM*„ Mk) = E 44/^fc), (B2) 



where the state index k now labels sites at which a 
nuclear spin has been raised or lowered, and with the 
help of the matrix elements: (I, m ± 1| I ± \I, m) = 



T m)(I ± m + 1). we have 



El 

3=1 



ml (mi 



±1) 



(B3) 



We assume the initial nuclear system is uniform, so that 
c± is independent of the site index fc, and for a large 
number of degenerate states g n 3> 1 that contribute to 
we replace the sum over weighting factors \ctj\ by an 
appropriate probability distribution. This gives c± = c± 
with c± defined in Eq. (|38[) of the main text. 



of I±(t) is not monotonic on the interval k = 1, . . . , N, 
where it has appreciable weight. We divide the sum 
into t subintervals of width Afc w y. The summand 
is then monotonic over each of the t subintervals, and 
the Euler-MacLauren formula gives a remainder R < jj 
when the sum over each subinterval is changed to an in- 
tegral. Adding the errors incurred for each subinterval, 
we find (for t > 1): 



I±(t) 



si le 



±it 



R(t). 



(C2) 



The remainder term \R(t)\ < so the corrections can 
become comparable to the amplitude of the integral itself 
when t « \J N/2. This represents a strict lower bound 
to the time scale where the continuum limit is valid for 
m = d = 2. 



Appendix D: PERTURBATION THEORY 



Appendix C: CONTINUUM LIMIT 

Here, we find a rigorous bound on corrections to the 
memory kernels, after we have changed sums to integrals. 
We consider the real-time version of the functions I±(s), 
given in (|37|l . with coupling constants for a Gaussian 
wave function in two dimensions (m = d = 2 in Eq. 
©): 



4N ^ 

k 



A 2 e ±iA k t/2^ 



-4, 



2e 



-k/N 



(CI) 



The Euler-MacLauren formula gives an upper bound to 
the corrections involved in the transformation of sums 
to integrals for a summand that is a smooth monotonic 
function of its argument. For times t > 1, the summand 



In this Appendix we apply standard perturbation the- 
ory to the problem of finding the electron spin dynamics. 
We do this to illustrate the connection between our per- 
turbativc expansion of the self-energy and the standard 
one, and to demonstrate the need for a non-perturbativc 
approach. 

We choose the initial state 



V // M0)IT) + e ? V 1 -'°T(0)U)) ®|n>, (Dl) 



where 



eigenstatc of h z and (5^)0 
(S+) = y/p^0)(l- Pl (0))e^. We 



|(Pt(0)-Pi(0)), w+/o 
then apply standard interaction picture perturbation the 

ory to evaluate (Sx) t , X = +,z. To lowest nontrivial 

(second) order in the perturbation TLy, we find 



J 



(S 



+/t 



at 



(S+) 



te 



iuj n t 



E ([*-]«* l h +]kn9k(t) + Ih+U [h-\ kn 9t{t)) 



1 

2 



(&>o + \ E K 1 - Mo)) [h + ] nk [h.] kn - pr(o) [h.] nk [h + ] kn ] I 1 " 



(D2) 

(D3) 
(D4) 



Ilk 



The expression for (S z ) t has been given previously;ii*i^ 
where it was noted that the perturbative expression for 
the transverse components (S+) t contains a term that 
grows unbounded in time (as above) . Inserting an initial 
nuclear state \n) with uniform polarization, performing 



the continuum limit, and expanding to leading order in 
— gives the final result, presented in Eqs. (|44fl and 145|l . 
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Appendix E: BRANCH CUT ASYMPTOTICS 
1. Long times 



grals for long times 



Here we give explicit expressions for the leading-order 
terms in asymptotic expansions of the branch cut inte- 



J 



KS(t ->oo) 



1 



(c+ +c„)7V 



[2(S 2 ) -A( c+ +c_)] — + 



in 
V 



N 



T (c++c-)[2(S z ) -A}-2pIA 



,±it 



Nc 



T I 



1 



ln 2 i 
1 



lni 



2 (S z ) T A ( — ± (c+ + c_) T c±21n2 



6' = 0, 
, b' ± 0, 
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Hn 2 t 



O 



(El) 
(E2) 
(E3) 



*tf(t->oo) 

iC£(t->oo) 



: ( C+ + c_)iv <5+>0 in7 + o (h77 
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2?re 



NCa 



(S+)o 
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iln 2 i 



1 



t\n 3 t 



b' = 0, 



iMc + + C-)-^{S + ) ±+o(±) 



(E4) 
(E5) 
(E6) 



2. High fields 

For asymptotically large magnetic fields, the in- 
dependence of the denominator term D(—x + j a ± irj) 
that appears in the branch cut integrals (Eq. i|74|H is 
dominated by the constant contribution ~ —iuj n , ex- 
cept at very large values of x, where it may be that 
D(—x + 7 Q ± irj) ~ 0. We expand the numerator 
N z (s^(x) ± irj) and denominator in retaining terms 
up to O (1) in the numerator and O ( in the denom- 
inator. Expanding to leading order in oc except 
where there is the possibility of a near-singular contribu- 
tion (D ~ 0), and assuming b' > 0, we find the branch 
cut integrals 



2niSC^x 



— 2+00 



dz 



ze 



z - z 



-2TTiS (Cf + Cf) x / dz 

Jo z — z o 



-2m5C 



x 



z+oo 



dzze zt 



(E7) 
(E8) 
(E9) 



where zq and xq are defined in Eqs. jHSjl and l|81(l . The 
coefficients C^are given by Eq. I|51|) . The sum over all 
three branch cut integrals can now be written in terms 
of two contour integrals 



a=(0,+,-) 



-27r^C5x / dz — 

J c „ z - z 



2m5C^ I dzze~ zt . (E10) 



C" runs clockwise from the origin to z = oo along the 
real axis, then returns to z = —i, enclosing the pole at 
z = zq. C" runs from z = i to z = i + oo, then returns 
along the real axis to z = 0. These integrals can be 
evaluated immediately by closing the contours along the 
imaginary axis. The sum of the contributions along the 
imaginary axis and from the residue of the pole at z = zq 
gives the result in Eq. I|79|) . 
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r sw is, e.g., the time taken to inject an electron into a 
quantum dot. 

It may be possible to measure the Overhauser field di- 
rectly by locating the position of the electron spin reso- 
nance (ESR) line, where the magnetic field compensates 
the nuclear Overhauser field. We have confirmed by exact 
diagonalizations on small (iV tot = 15) spin systems that 
the resonance is indeed centered at a magnetic field cor- 
responding to the negative nuclear Overhauser field, even 
for a nuclear spin system with p < 1. Alternatively, a state 
where all nuclear spins are aligned along the magnetic field 
can be generated by allowing the nuclear spins to relax in 
the presence of the nuclear spin-lattice interaction. 
Abe et all— have recently measured a pure Gaussian de- 
cay of the Hahn spin echo envelope with time scale given 
by the dipolar correlation time Tdd ~ 10~ 4 s for elec- 
trons trapped at phosphorus donors in isotopically en- 
riched 29 Si : P, where all silicon nuclei carry spin I = i. In 
contrast to the CPMG echo sequence, only a single mea- 
surement (a single echo) is made following each prepara- 
tion in the Hahn technique. We assume the echo envelope 
is the product of a Gaussian with time scale Tdd and a part 
/(2r) = 1 — 0(<5), 2r > tm, that gives the decay due to the 



contact hyperfine interaction: exp 



/(2r) 



O(S), for times tm < 2t <g Tdd- When 



_ l / _2r_ 

2 \T dd/ 

\ 2 

J <Si S, the dominant contribution comes from /(2t) 
at each echo of the CPMG sequence. 



